The von Neumann- Wigner theorem in quantum dot molecules 
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We show that electrons in coupled quantum dots characterized by high aspect ratios undergo 
abrupt density rotations when the dots are biased into an asymmetric confinement configuration. 
Density rotations occur with electron transfer to a single dot, and give rise to sharp variations of 
the exchange coupling between electrons as a function of inter-dot detuning. Our analysis based on 
exact diagonalization technique indicates that this unusual behavior is in agreement with the von 
Neumann- Wigner theorem that dictates the variations of the energy spectrum from the symmetries 
of the molecular states during the detuning process. It is also shown that the overall effect is 
quenched by the presence of magnetic fields, which by adding angular momentum to the electron 
motion affects the spatial symmetry of the molecular states. 

PACS numbers: 73.21. La, 73.21.-b 



The experimental observation of atomic-like proper- 
ties in manmade semiconductor quantum dots (QD) 
with feature sizes of a few hundreds to thousands of 
Angstroms has stimulated intense research in the last 
decade 0, 0, H, 0, [E II] • The existence of shell structures 
and the demonstration of Hund's first rule with shell fill- 
ing controlled by electric gates in QDs have provided the 
conditions for investigating tunable many-body effects in 
solid state systems 0, HI- As natural extensions of the 
concept of artificial atoms, coupled QDs resemble natu- 
ral molecules and offer enhanced many-body tunability 
by electrically varying the coupling barrier or distance be- 
tween the dots [9|, |l0( , and detuning the energy spectra 
of the individual QDs by appropriate gate biasing [ill ], 
which arc fundamental technical ingredients for realizing 
a quantum logic gate 12, 13| 

In this context, the symmetry of the electron wave 
functions is of crucial importance for the evolution of the 
energy spectrum of molecular systems under the influ- 
ence of perturbations. According to the von Neumann- 
Wigner theorem 14 1, energy levels cross for states that 
bear different symmetries, while they anti-cross for states 
with the same symmetry. This fundamental theorem has 
found successful applications in the spectroscopy of al- 
kaline salts [ll|, the interpretation of the Zeeman spec- 
trum of hydrogen molecules flfl ] , and more recently in the 
explanation of the magnetic anisotropy in thin metallic 
films 



17] 



Spatial symmetry effects can be mostly evidenced in 
double QDs with large aspect ratios, i.e., for which the 
confinement strength in the coupling direction is com- 
parable to or larger than that in the perpendicular di- 
rection. Hence inter-dot detuning is expected to induce 
significant conformational changes in the electron wave 
functions, which rearranges the electronic spectra of the 
individual QDs with drastic effects on the many-body 



interaction. 

In this Letter, we show that two electrons in coupled 
QDs with non-zero ellipticity experience density rotation 
and abrupt change of the ordering of their energy lev- 
els with inter-dot detuning. The charge density rotation 
occurs at specific detuning values, and leads to a discon- 
tinuity in the derivative of the exchange energy J with 
respect to the detuning bias, which is in agreement with 
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FIG. 1: (color online) Top panels: potential contour plots of 
coupled QD with r — 1 (left two) and r = 3 (right two). For 
each r, the panel on the left (right) is for Vl = Vr — 21 (Vl = 
31, Vr = 21) meV. (a) r = 1. Main panel: two lowest singlet 
(red, solid lines) and triplet (blue, dashed lines) energy levels 
as a function of detuning. U (V) indicates anti-crossing points 
for the two lowest singlet (triplet) levels. Lower inset: zoom- in 
region near the anti-crossing point of the singlet energy levels, 
(b) r = 3. Main panel and lowest inset same as in (a) but 
for r = 3. U' (V) indicates the anti-crossing (crossing) point 
of the two lowest singlet (triplet) energy levels Upper 
inset: three lowest single-particle energy levels in asymmetric 
confinement. In (a) and (b), the singlet energy levels are 
lowered by 3 meV for clarity. 
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FIG. 2: (color online) (a) Two-dimensional density plots for 
r = 1 at zero magnetic field. Columns (I), (II) and (III) 
correspond to Vl = 21, 29 and 42 meV, respectively. First 
(second) row, labeled S (T), is for the lowest singlet (triplet) 
states, (b) Same as (a) but for r = 3. Columns (I), (II) and 
(III) correspond to Vl = 25, 25.45 and 25.47 meV, respec- 
tively. Coordinates are shown in lower left panel and are the 
same for all panels. For r = 3, the peak separation along 
the x (y) direction in the density of the triplet state is ~ 60 
(~ 40) nm for column (II) [(III)]. 



the von Neumann- Wigner theorem [14|. We also show 
that the presence of magnetic fields enhances electron lo- 
calization in the QDs and breaks the spatial symmetry 
in the two-electron states, which quenches the density 
rotation effect. 

We model the coupled QDs with the ad hoc potential 
V(r) = -V L exp[-(x + d/2) 2 /R 2 x -y 2 /R 2 y ]-V R exp[-{x~ 
d/2) 2 /R 2 — y 2 /R 2 ], where d defines the inter-dot sepa- 
ration, R x and R y are the extensions of the Gaussian 
potential for each dot in the coupling direction x and 
perpendicular direction y. In this work, we take d = 60 
nm, R x — 30 nm, and define the QD aspect ratio as 
r = R y /R x . The difference between Vl and Vr gives the 
inter-dot detuning € = Vl — Vr. We fix Vr at 21 meV 
and increase Vl from 21 meV to higher values, which is 
equivalent to increasing e. Figure [TJ top panels show the 
potential contour plots for circular (r=l, left two panels) 
and elliptical (r=3, right two panels) QDs under symmet- 
ric and asymmetric biases. As e increases the coupled QD 
system becomes effectively a single QD centered at the 
lowest point of the confinement potential (the left dot 
center for both cases). 

We use numerical exact diagonalization of the two- 
electron Hamiltonian to solve for the system energies 
within the effective mass approximation. Details for the 
method and the material parameters are published else- 
where [18]. The convergence of the energies is tested by 



comparing the results using 8 and 9 harmonic oscillator 
basis states in each Cartesian direction (we use 9 states 
in each direction for this work). We ignore the Zeeman 
effect which simply lowers the triplet energy and J by 
~ 25 LieV/T at a given magnetic field. 

In Fig. [TJ bottom panels, we plot the two lowest singlet 
and the two lowest triplet energy levels for the coupled 
QDs with r = 1 [Fig. [Ha)] and r = 3 [Fig. [1(b)] as 
a function of detuning e (from here on, if not specified, 
"singlet" and "triplet" refer to the singlet and triplet 
states of lowest energy, respectively). All energy levels 
decrease as Vl increases because the single-particle ener- 
gies that give the main contribution to the two-electron 
energy are lowered with the electric potential. For both 
aspect ratios, r = 1 and 3, and at small e [21 < Vl < 29 
meV in Fig. [1(a) and 21 < V L < 25 meV Fig. Hps)], we 
observe singlet and triplet energies decrease with nearly 
the same slope as e increases. Then, the singlet energy 
slope becomes steeper as both electrons localize into one 
QD, which are indicated by points U and U' in Fig. [JJ 
For r = 1 and Vl > 34 meV, the slope of the triplet 
energy increases in value (yet smaller than that of sin- 
glet), while for r = 3 the triplet energy experiences a 
sudden change of slope at Vl — 25.46 meV [point V 
in Fig. UJb)], and with increasing Vl values afterwards, 
the spacing between singlet and triplet energies first de- 
creases, and then gradually increases for Vl > 27.2 meV 
(cf. J curve in Fig. [3]). 

The detuning effects on the singlet and triplet states 
can be further illustrated by inspecting the electron den- 
sity variations. In Fig. [2(a) we see that for coupled cir- 
cular QDs (r = 1) with increasing e electrons gradually 
move into the lower (left) QD for both singlet and triplet 
states, albeit more quickly for the singlet state. When 
both electrons are in the left dot [Fig. [2(a), column (III)], 
the singlet density has a single maximum corresponding 
to the state with nearly zero angular momentum [20j ]. 
while the triplet state exhibits two maxima along the QD 
coupling direction (x) as the constituent QDs are slightly 
more extended in the x-direction than in the y-direction. 
For r — 3, as e increases, the singlet density gradually 
localizes into the left dot with the electron density show- 
ing two maxima in the y-direction once Vl > 25 meV 
because of the relaxed confinement in that direction [cf. 
Fig. [2(b) columns (I) and (II)]. For the triplet state, how- 
ever, it is seen that the electron density abruptly changes 
from being spread over the two QDs with a higher peak in 
the left dot [V L = 25.45 meV, column (II) in Fig. [2(b)] 
to occupying only the left dot with two peaks of equal 
height in the y-direction [Vl = 25.47 meV, column (III) 
in Fig. [2(b)]. This abrupt transition (or "rotation") of 
the electron density occurs at Vl = 25.46 meV [point V 
in Fig. |2(b)]. 

In order to understand the drastic difference in the 
dependences of the energy levels on e between the 
QD configurations with r = 1 and r = 3, we com- 
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FIG. 3: (color online) Exchange coupling J as a function of 
Vl (Vr fixed at 21 meV) for QD aspect ratios r = 1 (black, 
solid), r = 1.5 (blue, dashed), r = 3 (green, dashed-dotted) 
and r = 5 (red, dotted) at zero magnetic field. The scaling 
factor for each curve is shown in parenthesis. 



pute the expectation values of the parity operator 
P) = {ty(xi,yi,X2,y2)\ty(-Xi,-yi,-X2,-y2)), and of 



the parity operator with respect to the y-axis (^P y 

{^(xi,yi,X2,y2)\^(xi,-yi,X2,-y2))- Here, * is the 
two-electron wavefunction, and {x\, y±), (X2, 2/2) are the 
coordinates of the first and second electrons, respectively. 
In the investigated detuning range, we find that for the 
singlet state, < P < 1 and P y = 1 for all r; while for 
the triplet state at r = 1, P < and P y = 1, as expected 
from the symmetry of the Hamiltonian. However, in the 
triplet state for r = 3 at Vl = 25.46 meV, our calculation 
shows that the value of sharply increases from —0.98 



to —0.02, while \Py/ changes sign from 1 to —1. This 
indicates that the parity of the singlet state with respect 
to the y-axis remains even, while the triplet wavefunc- 
tion parity changes from even to odd, thereby changing 
its symmetry along the y-direction. In fact, we find that 
as the inter-dot detuning is increased, the two lowest sin- 
glet states at any r have the same symmetry, and their 
energy levels anticross; the same is true for the two lowest 
triplet states at r — 1 (see Fig. [I}, whereas the energy 
levels of the two lowest triplet states at r > 1 [Fig. [ljb)] 
cross due to their different symmetries in agreement with 
the von Neumann- Wigner theorem [3]. Near the cross- 
ing point of the two lowest triplet energy levels [point V' 
in Fig. [Hb)], the analysis of the spectral function [2l| re- 
veals that the first excited single-particle state localized 
in the left QD has the p y -character [see schematic energy 
diagram in the upper inset of Fig. [ljb)]. Before the level 
crossing, this p y -\ike orbital is unoccupied so that the 
electronic states form a sp^-pair consisting of the s-like 
state mostly localized in the left QD, and the lowest anti- 
symmetric p^-like state spreading over both QDs. After 



the level crossing, it becomes energetically favorable for 
the triplet electrons to be in the same QD as an sp^-pair, 
even though the expectation value of the Coulomb inter- 
action in the triplet state increases from 2.0 meV before 
the transition sp x — > sp y to 2.9 meV afterwards. 

Figure displays the exchange coupling J as a func- 
tion of Vl, or effectively as a function of e = Vl — Vr. 
Each curve in Fig. [3] depicts the variation of J as the cou- 
pled QDs undergo the transition from the (1, 1) to (2, 0) 
charge states (the right end of each curve is chosen at a 
bias point where e is large enough to localize both elec- 
trons into the left dot). For small departure from equal 
biases (e > 0), as Vl increases from 21 meV, J main- 
tains a small and approximately constant value before it 
increases superlinearly, as was recently observed both ex- 
perimentally [H} and reproduced theoretically 22 1. For 
r = 1, J continues to increase monotonically with e as the 
two electrons strongly localize in one QD. In this case, the 
localization of the first excited single-particle state that 
has ^-character in the left dot occurs at larger detuning 
value than the ground state [points U and V in Fig. Q^a)], 
and J can only increase as the inter-level separation in 
the individual QD is larger than that in the coupled QD 
system. However, with increasing aspect ratios (r > 1.5), 
dramatic changes in the J behavior appear: a kink in J 
occurs at Vl = 30.70 meV with a sudden slope change for 
r = 1.5, while sharp maxima in J arise at 25.46 and 23.30 
meV for r — 3 and r = 5, respectively. In the latter cases, 
there is also a local J minimum at Vl — 27.20 and 26.20 
meV, respectively. For increasing r, electron localization 
in the left dot occurs at progressively smaller detuning 
(Fig. [3]), as it becomes easier for the system to lower 
its energy in a configuration where the two electrons are 
located at the opposite "ends" of the elliptical QD [Fig. 
H{b)] to minimize their Coulomb interaction. As a result, 
both singlet and triplet electrons move into the left dot 
at close Vl values [e.g., points U' and V' in Fig. \V[b) for 
r = 3] . Depending upon the difference between the biases 
required to localize singlet and triplet states into a single 
dot, the crossing between the two lowest triplet states 
associated with the electron density rotation is observed 
as a kink (r > 1) or a sharp maximum (r > 1) in the 
exchange coupling at intermediate values of detuning. 

Figure [4j main panel shows the ./-dependence on Vl 
for r = 3 in the detuning range corresponding to the 
(1, 1) (2,0) transition (24 < V L < 26 meV; V R = 21 
meV) for different magnetic fields normal to the xy-plane. 
We observe that the sharp discontinuities in dJ / dVL due 
to abrupt electron density rotation with detuning dis- 
appears for < B < 2 T, and evolves into a smooth 
J-increase with Vl- By adding angular momentum to 
the electron motion, the magnetic field mixes the p x and 
p y state parities, which affect the symmetry of the two- 
electron wavefunction, and quenches the abruptness of 
the density rotation. This symmetry breaking (from p y - 
symmetry to p x ~Py symmetry mixing) removes the cross- 
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FIG. 4: (color online) Main panel: J as a function of Vl for 
r — 3 at different magnetic fields. The magnetic field value 
is given on top of each curve. The curve for B — 0.1 T is 
indicated by an arrow. For clarity, curves at different B fields 
are shifted vertically in multiples of 0.03 meV. The dashed 
(dotted) gray line is a guide for the eyes tracing the localiza- 
tion point of the singlet (triplet) state at different magnetic 
fields. Inset: two lowest triplet energy levels at B = T (red, 
dashed) and B = 0.1 T (blue, solid). 



ing conditions between the two lowest triplet states to 
produce level anti-crossing, in agreement with the von 
Neumann- Wigner theorem (Fig. [4] inset). Our analy- 
sis also indicates smooth variations of the electron den- 
sity with detuning at small non-zero ( < B < 2 T) 
magnetic fields: for small e, the triplet electron density 
has a higher peak in the left dot than in the right dot, 
which, with increasing e, evolves into two peaks along 
the y-direction, while the peak in the right dot gradu- 
ally disappears. At higher magnetic fields (B = 3 and 
6 T), the localization of both singlet and triplet states 
in the left dot occurs abruptly with increasing detuning 
because of the strong magnetic compression of the elec- 
tron orbitals. We observe that for increasing magnetic 
fields the triplet state localization requires progressively 
smaller detuning [shown by the dotted line on Fig. SJb)], 
while the singlet state localization occurs at almost the 
same detuning [shown by the dashed line on Fig. EJb)] 
for all investigated magnetic fields [23j ]. 

Experimentally, the predicted variations in the two- 
electron density and exchange coupling can be obtained 
from the conductance responses of quantum point con- 
tacts adjacent to the dots, which are sensitive enough to 
register very slight changes in the Coulomb potential dis- 



13, 11 1 . One can 



tributions in the coupled QD system 
also detect the existence of the maximum in the exchange 
coupling from the coherent Rabi oscillations in the two- 
electron system by properly biasing the QD system 11]. 
Finally, electron imaging spectroscopy can be utilized to 
directly probe the electron density in QDs [24 1 . 



In conclusion, we show that inter-dot detuning in cou- 
pled elliptical quantum dots induces abrupt density rota- 
tions. We find that at zero magnetic field as the detuning 
increases, energy level crossing occurs for the two lowest 
triplet states bearing different (even/odd) symmetry. De- 
pending on the dot aspect ratio r, such a crossing results 
in a kink (r > 1) or sharp local maximum (r > 1) in 
the dependence of the exchange coupling on the detun- 
ing. For small nonzero magnetic field, the triplet energy 
level crossing in coupled dots with large aspect ratios 
evolves into anti-crossing due to the mixing of single- 
particle states bearing different spatial symmetry. For 
even larger magnetic fields, the localization of electrons 
by detuning is more abrupt due to the strong compres- 
sion of electronic orbitals. 
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